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Abstract 

Extensive variational computations are reported for the ground 
state energy of the non-relativistic two-electron atom. Several dif- 
ferent sets of basis functions were systematically explored, starting 
with the original scheme of Hylleraas. The most rapid convergence 
is found with a combination of negative powers and a logarithm of 
the coordinate s = r\ + r?,- At N=3091 terms we pass the previous 
best calculation (Korobov's 25 decimal accuracy with N=5200 terms) 
and we stop at N=10257 with E = -2.90372 43770 34119 59831 11592 
45194 40444 . . . 

Previous mathematical analysis sought to link the convergence rate 
of such calculations to specific analytic properties of the functions in- 
volved. The application of that theory to this new experimental data 
leaves a rather frustrating situation, where we seem able to do little 
more than invoke vague concepts, such as "flexibility." We conclude 
that theoretical understanding here lags well behind the power of avail- 
able computing machinery. 



* E-mail : schwart z@physics . berkeley. edu 



1 Introduction 



For thousands of years mathematicians have enjoyed competing with one an- 
other to compute ever more digits of the number tt. Among modern physi- 
cists, a close analogy is computation of the ground state energy of the Helium 
atom, begun 75 years ago by E. A. Hylleraas |J. 

Many have contributed incremental steps in this endeavor, flexing their 
computational muscle and ingenuity, often trying to use mathematical insight 
for advantage. The strongest line of theory has been to focus on the analytic 
properties of the wavefunction, especially following the 1954 work of V. A. 
Fock |2j which showed a weak logarithmic singularity at the three-particle 
coalescence. 

The recent work of V. I. Korobov || stands out for its simplicity and 
its success. His trial functions use the three scalar coordinates packaged as 
e -an-/3r 2 -7ri2^ w jth many sets of the complex nonlinear parameters cn,/3, 7 
selected in a quasirandom manner within specified intervals A, < «j < B ir , 
etc. With a linear combination of N=2200 terms of this type, grouped into 
four sets [Ai, Bj\, Korobov surpassed the best previous work by three decimals 
of accuracy in the Helium ground state energy; and more recently he went 
on to N=5200 and added four more decimal places. 

What struck me as surprising in Korobov's work was the fact that it 
seemed to ignore that earlier "wisdom" about analytic properties of wave- 
functions. His basis functions are, analytically, no different from the original 
Hylleraas basis; yet his computational results appear to converge so much 
more rapidly. This perception motivated the investigations reported below. 

Section 2 presents the new experimental data - systematic variational 
calculations using a variety of different basis functions that have been sug- 
gested over the years: these include negative powers, fractional powers and 
logarithms of the coordinates. Vastly different rates of convergence are ob- 
served, as illustrated in Figure 1. 

Section 3 offers a qualitative discussion and attempts to interpret this 
wealth of new data. Section 4 is a review, and an attempt to apply, the 
theoretical approach for understanding, at least semi- quantitatively, the ob- 
served convergence rates. We conclude that this state of the theory is far 
from satisfactory. 
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2 Experiments - Data 



Several different sets of basis functions were used in the standard variational 
calculations for the ground state energy of the Hamiltonian (with Z=2), 

H = -±[ V ? + V2] " Z/n - Z/r 2 + l/r 12 (2.1) 

and they are detailed below in terms of the Hylleraas coordinates, 

s = r 1 +r 2 , t = r 1 -r 2 , u = r 12 = \xi — x 2 \. (2.2) 

Table I presents summary results for the primary bases studied. The Ac- 
curacy of any energy value E is defined as Log[E*/ (E*-E)] where E* is our 
best estimate of the exact value. Thus, the numerical value of Accuracy tells 
how many correct decimal places there are in the calculated result. Gener- 
ally, we see that convergence was found to be more rapid as one progressed 
through this series, A,B,C,D,E,F. More detailed discussion of the results will 
be deferred to the following sections. 

• Basis A is just the original Hylleraas set: 

^ = Y.Ci, m , n e- ksl2 s l u m t n (2.3a) 
l,m = 0,1,2,3, n = 0,2,4,6, ... (2.3b) 

and we use the order parameter u = I + m + n to group the terms. We also 
designate a calculation of order uj to mean a basis set including all terms 
with / + m + n < uj. The total number of terms, N, grows asymptotically 
as a; 3 /12. This same organizational scheme is used for all the experiments 
listed below. 

• Basis B allows negative powers of s, which were introduced by Kinoshita 

§: 

e~ ks/2 s l {u/s) m {t/s) n - (2.4) 
and we can rearrange the terms as 

e~ ks ' 2 s l (u/s) m+n P n (t/u), (2.5) 

using the Legendre polynomials to get the selection rule An = 0, ±2. This 
technique follows the work of Goldman || and leads to more efficient use of 
computer space and time. 
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• Basis C allows fractional powers of s, which were first introduced by 
H. M. Schwartz §: 

(M 1/2 ) e~ ks/2 s l u m t n (2.6) 

which doubles the size N of the basis for each order uj. (The s 1 / 2 term is 
omitted for uj = 0.) 

• Basis D uses the logarithm of s, first introduced by Frankowski and 
Pekeris [0: 

(ljn(s)) e~ ks/2 s l u m t n (2.7) 

(The ln(s) term is omitted for u = 0,1. The values of N shown in Tables I 
and II are two units off for bases D and F.) 

• Basis E has both negative powers and fractional powers of s. 

(1, s 1/2 ) e~ ks/2 s l (u/s) m (t/s) n . (2.8) 

• Basis F has both the logarithm and negative powers of s. 

(1, ln{s)) e~ ks/2 s l {u/s) m (t/s) n . (2.9) 

For each of the basis sets described above, the scale parameter k was 
varied somewhat to find the lowest energy at each order. For the Table I 
data on Basis A, k varied from 5.0 to 8.2; for Basis B, from 3.8 to 5.9; for 
Basis C, from 4.7 to 6.6; and for Basis D, from 4.6 to 6.6. For Bases E and F, 
it was found that the optimum value of k stayed close to 2.0 for the mid-size 
and larger orders, so k was fixed at this value for all the data shown. 

In two cases I replaced the set of functions e~ ks ^ 2 s l by the set e~ aiS , using 
Korobov's quasi-random method for selecting the (real) nonlinear parameters 
a in a single group. The results for these experiments are shown as Bases B' 
and E' in Table I. 

Several variants of these basis sets were also explored briefly but discarded 
when they appeared less effective, as functions of N, than their counterparts 



above. Among these were the use of: 

Negative powers of s and u {u/s) m (t/u) n (2.10a) 

Third roots of s instead of the square root (2.10b) 

Two or more powers of ln(s) (2.10c) 

The coordinate r = Jr 2 + r 2 (2.10d) 



The coordinate R — \x\+xi\ (See the Appendix.) (2.10e) 
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Table I. Energy Accuracy - Summary data 



UJ 


N 


A- Basis 


B-Basis 


B' -Basis 


9 


125 


7.9 


8.7 


9.4 


12 


252 


8.7 


10.2 


10.7 


16 


525 


10.4 


11.7 


12.3 


21 


1078 


11.6 


13.2 


13.9 


27 


2135 


12.8 






UJ 


N 


C-Basis 


D-Basis 




7 


139 


9.5 


9.6 




9 


249 


11.2 


11.3 




12 


503 


13.4 


13.5 




16 


1049 


15.8 


16.0 




21 


2155 


18.3 


18.5 




UJ 


N 


E-Basis 


E' -Basis 


F-Basis 


7 


139 


9.5 


9.8 


9.4 


9 


249 


11.7 


11.6 


11.5 


12 


503 


14.5 


14.5 


14.3 


16 


1049 


18.2 


18.3 


18.1 


21 


2155 


22.5 




22.4 


27 


4269 


27.2 




27.6 


34 


8093 


29.7 




33.0 



Accuracy = Number of correct decimals 



In Table II are the detailed results (in a compact format) for the two 
bases - E and F - that show the most rapid convergence. One quantitative 
measure of the rate of convergence is the Ratio of successive differences and 
this is also shown in the table. 

Technical notes. For these computations, I wrote a set of subroutines for 
multiple-precision arithmetic (in C), eschewing more professional packages 
which are available. The last row of data in Table II used 101 decimals of 
precision and took one week running on an otherwise idle desktop computer 
equipped with a 300 MHz processor and 320 MB of memory. 
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Table II. Calculated results with the two best bases 



CO 


N 


R-Rasis Rnprcnps 


R /}t,ins 

X LCI 11UO 


R-Rasis RnprPHPs 


R f\ tins 

X LCI U Ivu 


A 


43 


-2 90372 1 83 

Zj.<JVJU I Zj lUa 




-2 90371 941 

Al.L/WtJ 1 X J11 




U 


67 


-2 90372 42300 


1 7 9 

X 1 . U 


-2 90372 41 5 


22 5 

Z, Z, . ' J 


6 


QQ 


" 4S64S 




-2 00^72 4^61 


14 2 

It: . Zj 


7 




43762 2 

t: r .J | VJZj 


15 8 


" 43758 7 


14 


8 


1 8Q 


43769 759 

t:-J ( VJ.J 1 UtV 


14 3 


43769 382 

( VJ.J tJUZj 


12 2 

X Zj . Zj 


Q 


24Q 


43770 28348 


10 5 


43770 25283 


1 1 2 

X X . Zj 


10 


321 


43770 33352 


7 32 


43770 33068 


8 35 


11 


405 


" 34036 


9 37 


" 34000 4 

■J 1 VJ VJ VJ 1 


8 80 


1 2 


503 


341 09 6 

tJTlUJ VJ 


8 30 


341 06 294 


8 94 


1 3 

X U 


61 5 

VJ X U 


S41 1 8 444 

ijtIIU XXX 


8 63 

O . VJU 


341 18 1 3633 


9 02 

U . WZj 


14 

X t: 


743 


S41 1 9 46685 

(Jt:11 J T:VJVJ(JtJ 


8 79 


S41 1 9 44847 

tjUlJ 1 IUt: 1 


10 1 

1U.1 


1 5 

X U 


887 


341 1 9 58229 


8 06 


341 19 57846 


7 33 


16 


1040 


" 59667 


10 4 


" 59620 

XJU VJZj VJ 


10 


17 

X 1 


1229 


59806 


6 32 

V . - J Zj 


59797 


5 90 

?J . l/U 


1 8 

X o 


1429 

XtiZj ,J 


59828 5 

U<J(JZj(J tj 


10 5 


59827 3 


9 38 

U . Uu 


1 Q 

X <v 


1 649 

XVJt:<J 


59830 654 


5 12 

tj . X Zj 


59830 456 


5 50 


2D 


1 891 

X O U X 


59831 06419 

u J(J«J1 VJVJt:X<J 


9 69 

<J . VJ.J 


59831 03831 


8 89 


21 

Z, X 


21 55 

Zj X xjxj 


59831 10650 

tj JU«J1 XVJVJUVJ 


5 08 

*J . VJ(J 


59831 10381 


6 18 

VJ. X(J 


22 


Zj X 1U 


" 1 1482 


9 14 


" 1 1 442 

X X 1 XZj 


8 22 


23 

ZjU 


2755 

Zj 1 UtJ 


1 1573 8 

X XtJ I u o 


5 62 

\J . VJZj 


11571 6 

lltj 1 1 VJ 


7 25 

1 . ZjU 


24 

Zjt: 


3093 


1 1 589 997 


7 94 


1 1 589 408 


6 66 

vj . vjvj 


25 

ZjU 


3457 


1 1 592 03902 

VJUiJVJZj 


6 18 

VJ. xo 


1 1592 08081 


8 87 


26 

Zj VJ 


3849 


1 1 592 36947 


5 29 

5j . ZjiJ 


1 1592 38154 

X X \J ■_./ Z. *J<JXtJt: 


4 79 


27 


4269 


" 43186 


4.80 


" 44444 


11.0 


28 


4719 


44484 


2.97 


45017 4 


3.57 


29 


5199 


44922 


2.88 


45177 752 


13.1 


30 


5711 


45074 35 


2.32 


45189 95689 


3.00 


31 


6255 


45139 97 


2.29 


45194 02040 


14.8 


32 


6833 


45168 66 


2.16 


45194 29566 


2.76 


33 


7445 


45181 95 


2.11 


" 39521 


15.1 


34 


8093 


45188 24 


2.06 


40182 


2.78 


35 


8777 


45191 29 


2.02 


40420 7 


13.5 


36 


9499 


45192 80 




40438 342 


3.11 


37 


10259 






40444 00495 





extrapolate: E* = -2.90372 43770 34119 59831 11592 45194 40444 6 
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3 Experiments - Discussion 



Figure 1 provides a visual comparison of the convergence rates for the dif- 
ferent variational basis sets (A - F), plotting Accuracy versus the Log of N, 
the number of basis functions used. I found a number of surprises in these 
results. 

Surprise 1. Basis B (negative powers of s) shows a significant improvement 
over basis A (the original Hylleraas basis). 

Surprise 2. Bases E and F do a great deal better than any of the others. 
Somehow, the benefits of B and C (or B and D) are cumulative. 

Surprise 3. The performance of C and D are nearly identical, as with E 
and F (until we reach very high orders.) 

Surprise 4. The performance of basis E drops off dramatically after 
uj = 27; but basis F keeps up its rapid convergence, although with marked 
oscillation, as seen from the Ratios in Table II. 

The surprising performance of Korobov's basis has already been noted: 
his published results are shown by the line labelled with the letter "K" in 
Figure 1. 

The fact that Basis C performs a lot better than Basis A was not a 
surprise, since earlier work || had already shown that. Ditto for Basis D M. 
Also, the smallness of the gain shown by basis B' over B (and by E' over 
E) is as expected, based upon the analytic equivalence of exponentials and 
power series. 

While I cannot explain the surprising results, I can readily offer sugges- 
tions on how one might interpret them. The virtue of Bases C and D lies in 
providing more flexibility to the "radial" behavior of the wavefunction (coor- 
dinate s); while that of the negative powers lies in providing more flexibility 
in the "angular" behavior (ratios u/s and t/s). The similarity between C and 
D (and between E and F) indicates that the precise analytic behavior in the 
"radial" coordinate is not important - any flexibility will do - until one gets 
to the very high orders. 

This appeal to "flexibility" is just armwaving; it lacks any mathematical 
foundation. Such appeal to flexibility is also the best way I know to un- 
derstand the success of Korobov's calculations: his work seems akin to the 
"adaptive" techniques used in numerical integration, where one puts addi- 
tional mesh points into any region that shows a slower rate of convergence. 

In varying the scale parameter k, I most always found a simple minimum 
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in the energy. However, in one case - Basis Eatc<j = 21-a more complex 
behavior emerged: see Table III. While the variation shown here is not very 
great, this does raise the general question of how effectively one may search 
for the minimum of a complicated function of many nonlinear parameters. 
This is a possible source of worry in using Korobov's technique, especially 
when it comes to error estimation and extrapolation. 



Table III. Double minimum in E(k) for Basis E at u = 21 



k 


1.7 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


2.4 


E 


10567 


10697 


10691 


10650 


10635 


10640 


10636 


10603 


These num 


3ers for the energy E fo 


low the first 20 decimal places. 



4 Theories 



The first lesson in analysis of atomic wavefunctions concerns the two-particle 
cusps: linear behavior as any one of the coordinates r±, T2, or r\2 goes to zero. 
All of the basis functions studied in this paper are correct in that regard; we 
are concerned here with what comes next. 

Take the Hylleraas expansion ( ^.3a| ) and put it into the Schrodinger equa- 
tion ( |2.1|) . Then collect the coefficients of each monomial in s, u, and t and 
set that equal to zero. Early in this infinite set of algebraic equations for the 
expansion coefficients C^ m ^ n one finds the following inconsistent equations 



Ci,o,o + ZC-ofifi — 
— 2Ci ) i i o + Ci,o,o = 
4Ci,i,o — Ci,o,o — 0. 



(4.1a) 
(4.1b) 
(4.1c) 



This contradiction in the Hylleraas basis was discovered by Bartlett, Gib- 
bons and Dunn || in 1935 and it led them to consider an alternative expan- 
sion: one that involved logarithms of the hyperradius r = Jrf + r\. Later, 
Fock |J independently developed a systematic expansion of the wavefunction 
with such terms. 



In 1962 this author |]TU[ developed a general theory about the conver- 
gence rate of variational calculations, based upon analogy with least-squares 
fitting of functions and one- dimensional model problems. This line of anal- 
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ysis was expanded by others [D] and in some cases given a more rigorous 
mathematical basis. | 12| 



That theoretical work led directly to the idea that the convergence rate 
in Hylleraas-type calculations for the Helium ground state was controlled by 
the Fock logarithmic singularity; and the semi-quantitative analysis seemed 
to fit the available data. It also led to the successful exploitation of the 
fractional power basis C |§. Shortly thereafter, Frankowski and Pekeris 
[0] took logarithmic terms explicitly into their trial functions and this also 
seemed to confirm the importance of the Fock behavior. 

However, when Kinoshita [|J considered the systematic use of negative 
powers - in the form ( [2.10a|) , not ( [2.4|) - he found that there was no contra- 
diction of the type noted above. I have confirmed that this conclusion holds 
also for the basis (|2.4|) used in the current work. 

Thus, one might be drawn to believe that the logarithmic singularity is not 
an intrinsic property of the He atom wavefunction itself but rather results 
from a particular analysis that relies on the six-dimensional coordinate r. 
This idea may be dispelled by reading the work of Myers et al ||13| , where 
they take one for a walk around the three-particle coalescence and show that 
the (finite) discontinuity in the local energy disappears when one includes 
the full set of terms that accompany the Fock logarithm. 

This approach helps us understand some other experimental results. We 
noted earlier that inclusion of negative powers t/u gave poorer results than 
t/s. Each of these ratios shows a (finite) discontinuity when one walks around 
the place where the denominator vanishes - something that the correct wave- 
function should not allow. In the case of (t/u) this error occurs along a line, 
while in the case of (t/s) this error occurs only at a point. A similar situa- 
tion appears in the work of Goldman ||, who holds the previous record for 
basis size (N=8066). His use of the coordinates r<,r> led to a very efficient 
computational scheme, but demonstrates relatively slow convergence. (See 
the point "G" in Figure 1.) This may be attributed to the discontinuity of 
his basis functions on the 2-dimensional surface ri = r 2 . 



4.1 Fitting the Data 

The goal of a good theoretical understanding should be the ability to pre- 
dict or to explain, at least semi-quantitatively, the observed rates of con- 
vergence for systematic variational calculations with different basis sets. In 
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my earlier work | ID| , the attempt to do this was based upon analogies with 
one-dimensional model problems, doing least-squares fit with appropriate 
orthogonal bases to represent functions with various types of singularities: 



Minimize 



n— 1 



p(x)dx [f(x) - CiUi(x)f 



i=0 



Ci = J p(x)dxf(x)ui(x) 
Error « (C n ) 2 . 



(4.2a) 

(4.2b) 
(4.2c) 



For one example we find: 

f(x) = x u Inx, 
C n 

and an alternative example is: 



p{x) 



X 



on the interval (0,1) 



1/n 1 



(4.3) 



/(*) 



x' 



Inx, 



p(x) 



x^e 



on the interval (0,oo) 



(4.4) 



The difference in convergence rates for these two examples may be understood 
qualitatively as follows. The basis functions x n e~ x peak at x = n. Therefore, 
at higher n these basis functions on the interval (0, oo) get farther and farther 
away from the singularity, which is at x = 0. One may improve the situation 
by using basis functions x n e~ kx , where k is a scale parameter that may grow 
as one proceeds to higher orders. I do not have a quantitative theory for this 
result but it is qualitatively relevant to the current study. 

In my 1962 work I applied this simple modeling to the He atom problem, 
identifying the Fock term r 2 lnr as the dominant singularity which is neglected 
in conventional Hylleraas coordinates. This led me to predict a convergence 
rate formula, 

E(lu) - E(lu - 1) ~ const. fu? (4.5) 

and I estimated that p should be between 5.5 and 10, due to uncertainties 
in replacing the real 3-dimensional problem with the one-dimensional model. 
The then best results with Hylleraas variables (work of Pekeris, 1959, up 
to order 21, using a cleverly orthogonalized basis) fit the convergence rate 



formula (|4.5|) with a value of p between 7 and 8. This was good confirmation 
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of the theory. The extended computations reported here (Basis A data in 
Table I) fit the convergence rate formula ( |4.5|) with a value of p which varies 
from 7, at the lower orders, to a value about 12 at the higher orders. This 
improvement is probably due to my allowing the scale parameter k to vary, 
which was not done in the earlier work. 

Also, in 1962, I introduced the half-powers of coordinate s, explicitly for 
the purpose of increasing the convergence rate, following this theory. That 
was successful, with the observed value of p approximately doubled to 14 or 
15 at uj < 8. The extended computations reported here (Basis C data in 
Table I) are fit to values of p which vary from about 16 to 21. Again, this 
is fairly good confirmation of the theory; and again we acknowledge some 
improvement by allowing the scale parameter to vary. 

Following that earlier theory one would certainly not expect Basis D to 
converge at the same rate as Basis C - but this is exactly the behavior we 
have found in the present experiments. 

What can I say about the observed convergence rate of Basis B, intro- 
ducing negative powers into the Hylleraas functions? The data in Table I are 
fit with a value of the exponent p around 13. I do not understand this but 
will only offer a guess that it may have to do with fitting the complex "angu- 



lar" behavior around the Fock singularity, which was described in Ref. [pT5 
Maybe this is connected with the difference in convergence rates noted above, 
in Eqs. ( f4.3| , jO]) , for the model problems on (0,1) and on (0,oo ). 

Finally, look at the results for Bases E and F. The data in Tables I and 
II are fit with values of the exponent p which grow from the 20's to the 40's 
in the middle range of uj; at the top end, the data for Basis E drop to around 
p=25, while the data for Basis F climb to about p = 65. I am at a loss to 
explain these large exponents following the former analysis. 

An alternative to the power law convergence rate formula ( [4.5[ ) is the 
exponential rate formula 

E(u) - E(uj - 1) ~ const. (a) w (4.6) 

which one could expect from a model fitting problem that involved no singu- 
larities at all. For example, expanding e~ ax in terms of x n e~ bx would yield 
the formula ( |4.6| ) with a = (f^pf) 2 . If one plots the data for Basis F (Log of 
increments in E vs. uj), it does look close to a straight line; and the smoothed 
data in Table II may be fitted with a value of a in the range 0.13 - 0.16 for 
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to > 16. If one looks at the asymptotic behavior of the He wavefunction as 
ri goes to infinity, the behavior in r2 should be as e~ Zr2 with Z = 2. The 
trial functions I used for this basis have the exponential envelope e - fc ( r i+ r 2)/ 2 
with k = 2. Using the formula quoted above, this model gives us the param- 
eter a as (fqrj-) 2 = 0.11. This looks like a fairly good fit to the data; but 
accepting this explanation would lead us to doubt the relevance of the Fock 
singularities for the He wavefunction. 

Also, I know of no published theoretical attempts to explain the excel- 
lent convergence found by Korobov with his highly nonlinear fitting of the 
trial wavefunction. John Morgan has suggested (in private communication) 
that Korobov's approach may be likened to the work of fitting the Hydrogen 
radial wavefunction with a set of gaussians, using "floating exponents" |T3 



This sounds plausible, but at present it is just more handwaving about "flex- 
ibility". 

I conclude that theoretical understanding of the convergence of variational 
calculations on the two-electron atom is far outstripped by the raw computing 
power of available machinery. 

Some may ask if any of this is really relevant to current issues in physics. 
One response is to point to high accuracy measurements performed on atomic 
systems which may check the current theories of fundamental particles and 



interactions. A recent paper by Pachucki and Sapirstein [JT5] aims to deter- 
mine the fine structure constant to a few parts-per-billion. This is based 
upon measurements of the 2 3 Pj states in Helium and detailed calculations 
that rely upon a Korobov-type representation of the atomic wavefunction. 
Then, again, all this may be nothing more than an expression of 7r-envy. 
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Appendix: Integrals 

Integrals of the following type were needed in the calculations reported 
here: 

dse~ s s p (ln(s)) q . (A.l) 
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There is a simple recursion formula on the index p; and for the minimum 
values of p I used a particular technique of numerical integration. (See Ref. 
IPI-) First change variables, s = exp(y); then use the simple rule, 

f(y)dy*h £ firth). (A.2) 

-°° n=-oo 

The summation is truncated when terms are smaller than the desired accu- 
racy; and the answer converges exponentially as the interval h is decreased. 

For the two-electron atom, one can evaluate the most conventional inte- 
grals from the formula, 

d 3 x ± fd 3 x 2 e" ari e~ br2 e" cri2 1 

' " (A.3) 



An J An n r 2 r i2 (a + b)(b + c)(c + a)' 
and derivatives of this simple result with respect to the parameters a,b,c. 

In exploring more complicated functions, I was able to find another simple 
formula for the following integral, which involves R— \x% + x 2 \. 

r d 3 x x r d 3 x 2 e~ ari e~ hT2 e~ cri2 e~ dR 
J An J Air r± r 2 r 12 R 

1 

(a 2 + b 2 - 

To derive this, insert the Laplacian operators into the middle of the integral 
and let them work both ways. It appears that one could almost deduce these 
results ( |A.3| ) and ( |A.4| ) purely by arguments of analyticity and symmetry. 
Consider, for example, how the integral behaves as r\ — > 00: by counting 
powers one sees the nature of the singularity as (a+c), or (a+c+d), goes to 
zero. 

As noted earlier, using this variable R in the He trial wavefunction did 
not produce good results - as one might expect since it introduces a spurious 
cusp when the two electrons are on opposite sides of the nucleus. I have, 
nevertheless, recorded the above information here in case it might be useful 
to others. 

The result ( |A.4|) can be generalized with R = \[ixi + vx 2 \- 



- (A.4) 

1 , (a + b + 2c)(a + b + 2d) 

/ 71 

- 2c 2 - 2d 2 ) 2(a + c + d)(b + c + d)' 
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